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Motivated by recent neutron scattering experiments on the cuprate superconductors, we present 
a phenomenological framework describing the dynamics of collective spin excitations coupled to 
charge/bond order fluctuations. Our quantum lattice model contains two order parameter fields, 
and can capture spin excitations both in broken-symmetry states with static lattice modulations, 
as well as in homogeneous states where the charge/bond order is fluctuating. We present results for 
different types of static charge/bond order, namely site- and bond-centered stripes, and plaquette 
modulation. 
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I. INTRODUCTION 

For certain cuprate high-temperature superconductors 
it has been estabhshed that incommensurate spin and 
charge correlations, commonly referred to as stripes, ap- 
pear over a significant range of the phase diagramiSii. 
The role of these stripes both for superconductivity itself 
as well as for various anomalies in the normal state has 
been discussed extensively^*^*SiL2i2ii£, but is at present 
not fully understood. 

Recent neutron scattering experimentsii^iSiiiii have 
mapped out spin excitations in various cuprates over a 
large range of energies. Tranquada et ali^ investigated 
the excitation branches in La2_a;Baa;Cu04 at a hole dop- 
ing of a; = |, which displays static spin and charge order, 
up to energies of 200 meV. They found the high-energy 
part to be well described by the spectrum of spin lad- 
ders, pointing towards bond order in the material^'^-^®. 
Remarkably, other cuprate families like YBa2Cu306+j; or 
Bi2Sr2CaCu208+5, where charge order remains dynami- 
cally fluctuating, show a very similar spin excitation spec- 
trum at elevated energies. (At low energies, the spectrum 
depends on whether the system shows static order; fur- 
thermore, the strength of electronic quasiparticle damp- 
ing is different between different compounds, limiting the 
wavevector range where sharp single-particle excitations 
can be seen.) These experiments raise the question of 
whether there exists a dynamic spin response which is 
universal among all cuprates: namely a high- intensity 
("resonance") peak at wavevector (tt, tt) and energy be- 
tween 20 and 50 meV, with both downward and upward 
dispersing branches of excitations. 

The purpose of this paper is to develop a unified 
theoretical framework for describing spin excitations in 
the presence of both static and dynamic charge order. 
We will do this using a phenomenological lattice model. 
Our approach differs from previous phenomenological 
theoriesiiiSiiS in two important ways: (i) These theo- 
ries took the continuum limit for spin modulations in the 
vicinity of the incommensurate spin-ordering wavevec- 



tor; such a theory is not expected to be valid near the 
resonance peak at (7r,7r). We will instead use a lattice 
model which is valid both at (tt, tt) and near the incom- 
mensurate wavevector, and for a wide range of energies. 
(ii) Our model explicitly selects colUnear spin correla- 
tions, in contrast to the previous approaches^^--'^^ where 
the distinction between coUinear and spiral spin correla- 
tions only appears in higher order terms which are not 
fully accounted for in existing calculations. 

Our model is an extension of recent microscopic spin- 
only models for stripe phasesi^ii^ which have quite suc- 
cessfully modelled some of the neutron scattering dataiii 
One of our purposes is to judge the significance of this 
agreement between theory and experiment: in particu- 
lar, we wish to delineate the range of models which are 
compatible with the data, including those which cannot 
be described by simple spin-only Hamiltonians. Further- 
more, our approach can be naturally extended to the case 
where the charge order is dynamic, rather t han th e static 
order needed to define the models of Refs. Il5ll6i 

We will model spin and charge fluctuations on a phe- 
nomenological basis using Landau order parameters for 
both. The coupling between the two orders can shift the 
minimum energy of the spin excitations from (tt, tt) to 
the incommensurate wavevector dictated by the charge 
order. On a microscopic scale, the influence of the 
charge order on the spin sector can be understood in 
terms of both spatially modulated spin densities and spa- 
tially modulated couplings, as found in the models of 
Refs. Il5ll6i Our Landau-like theory will be formulated 
directly on the underlying square lattice, which will al- 
low us to capture lattice effects such as the differences be- 
tween site-centered or bond-centered order. Importantly, 
the spin sector described by our theory will be strongly 
fluctuating, i.e., we are far from the semiclassical limit 
described by spin wavesi^*^'^^. Microscopically, these 
fluctuations can arise from the tendency to dimerization, 
i.e., from bond order, which is present in the undoped 
paramagnetic parent Mott insulator^'^^ . (Note that bond 
order can occur both in site-centered and bond-centered 
stripe states.) 
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II. QUANTUM LATTICE MODEL 

We assume a dominant antiferromagnetic interaction 
between the spins, and so model the quantum spin fluc- 
tuations by a standard vector tp'^ Landau theory for the 
antiferromagnetic order parameter at the commensurate 
wavevector Q = (tt, tt). So on a square lattice of sites, j, 
we parameterize the lattice spins by 



(1) 



where a = x,y, z. In the absence of any coupling to 
charge/bond order, we assume that the dominant spin 
fluctuations remain at the commensurate Q; we find be- 
low that this feature is important in obtaining a reso- 
nance peak at Q. The effective action for these commen- 
surate spin fluctuations has a familiar form: 



So = 



1 / dipja 

2\ dr 



(2) 



In principle, there should also be Berry phases in the 
quantum spin action, but we assume that they have av- 
eraged out to zero: this is expected to be valid in the 
compressible superconducting states, or in the incom- 
pressible Mott insulators with an even number of elec- 
trons per unit cell, but likely not at the quantum critical 
point between such phases. 

Now we include the effect of charge/bond order. This 
we represent by the complex continuum fields (/>a;^j,(r, t) 
which measure the amplitude of charge order at the 
wavevectors K^; = (7r/2,0) and Ky = (0,7r/2) - this is 
the dominant ordering wavevector of the Mott insulat- 
ing state at 1/8 doping, and is therefore the appropriate 
reference wavevector for our considerations: we will show 
below how deviations in the charging ordering wavevector 
from ii.x,y in the superconducting phases (or at non-zero 
temperatures) can be easily built into our formalism. Us- 
ing the complex order parameters (j)x,y, it is convenient 
to define the real field 



(3) 



and similarly for Qy. For r on the sites of the square 
lattice, the Qx,y are measures of the charge density mod- 
ulation on those sites. On the other hand, for r on the 
links of the square lattice, the Qx.y is a measure of the 
local bond order: this is determined by the modulation 
in the local pairing amplitude or exchange energy. With 
these physical interpretations at hand, we can write down 
the following couplings between the spin and charge fluc- 
tuations 

Sx^ dT^^XlQx{rj)(p'^j^+ X2Qx{rj+x/2)iPja^j+x.a 
'' 3 

+ X4Qx{'rj+y/2)fjafj+y,a\ (4) 



with four independent couplings constants Ai_4; the 
same couplings will appear in the corresponding Sy. No- 
tice that Ai implements the correlation between the on- 
site charge density and the amplitude of the spin fluc- 
tuations, while A2-4 ensure that the effective first- and 
second-neighbor exchange constants controlling the spin 
correlations modulate along with the bond order. 

If we now take (l3x,y — constant, then the action SQ-\-Sx.y 
represents our general theory for quantum spin fluctua- 
tions in a background of static charge/bond order. The 
site-centered case has (j)x — I, and the bond-centered case 
has (px = e^'^l^ . Two-dimensional (plaquette or checker- 
board) order will have both (px and 0y non-zero. At the 
Gaussian level (u = 0), the problem is quadratic in the 
Lp flelds, and can be solved by diagonalizing a matrix of 
size {N^Ny)'^, where N^^y describe the size of the unit 
cell (in the charge sector). 

For large enough A couplings, the minimum energy of 
the if fluctuations will be shifted away from (tt, tt), as ob- 
served in experiment. Notably, the restriction to real Lp 
implies that the spin order remains coUinear. For small 
mass s, the spin order can condense as usual, and fluc- 
tuations around the condensate will lead to low-energy 
Goldstone modes. 



A. Fluctuating charge order 

The above formalism is designed to allow easy exten- 
sion to the dynamic charge order case, which is likely rel- 
evant for YBa2Cu306+y and Bi2Sr2CaCu208+5. We use 
a continuum formulation to describe the charge fluctua- 
tions, with the following general action consistent with 
all the symmetries of the lattice: 



\dr(t>x\ +\dr(l)y\ +cl\dx(t>xY 



+ C2\dy4>x\ +C^\dy4>y\ 1 9^; | +i5(l)ldx4>x 

|2 I Ij, |2\ I „, /ij, |4 I ij, |4\ 



-lO0:.dy(t>y + Sl {\(t)x? + M) + Ul (10:, I* + \4>y\'^ 



(5) 



Note especially the term proportional to 5: it is generi- 
cally present, and when the strength of the charge order 
is weak it ensures that the dominant charge order fluc- 
tuations are at an incommensurate wavevector unequal 
to ^x,y However at low temperatures (T) in an insulat- 
ing state, the 'lock-in' term proportional to w eventually 
dominates, and selects a commensurate charge-ordered 
state with wavevectors equal to 'i^x,y The quartic v 
term determines whether the long-range order will be 
one-dimensional ( "stripe" ) or two-dimensional ( "plaque- 
tte", "checkerboard"). 

The effect of charge fluctuations on the spin excita- 
tion spectrum can now be determined by self-consistently 
computing the frequency and momentum dependence of 
the ifinSelf energy to second order in the A's, similar to 
Ref. [2^ A full calculation along this line is beyond the 
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scope of this paper, we expect that in the Hmit of small 
mass Si the results will be similar to the ones obtained 
in the static (j) theory. 



B. Coupling to phase fluctuations 

An advantage of the present phenomenological formal- 
ism is that is allows easy extension to include couplings to 
other collective modes. If the cuprate compound is a su- 
perconductor, then there is an additional mode associate 
the fluctuation of 9, the phase of the superconducting or- 
der. As in Ref. I23. the most relevant coupling of 6 to the 
spin fluctuations is 



80meV 



200 meV 



(6) 



The action for 9 fluctuations can be generally written 

„„25 



Se 



(Pkduj 



(7) 



where fc is a wavevector and u is an imaginary frequency, 
and the parameter a is determined by the nature of the 
Coulomb interaction: for screened short-range interac- 
tions fT = 0, while for in- plane 1/r Coulomb interactions 
with independent layers a — 1. In a paramagnetic state 
where excitation forms a sharp S = \ 'triplon' ex- 
citation, the coupling to phase fluctuations will induce 
damping in the triplon spectral function. At T = 0, and 
at the bottom of the triplon band, it is not difficult to 
compute from Eqs. Ht)l7|l that the imaginary part of the 
triplon self energy is ^ (e — A)^('*~'^'/(^~'^) at the bottom 
of the band (e > A is real frequency and A is the spin 
gap): this arises from the 'radiation' of 9 excitations by 
the triplon. 



III. RESULTS FOR STATIC CHARGE ORDER 

Let us now present a few result for the case of static 
charge/bond order, calculated in the Gaussian approxi- 
mation of So + Sx + Sy. We directly calculate the T = 
susceptibility as measured in inelastic neutron scatter- 
ing. In the following, we restrict our attention to the 
one-particle contributions, as the multiparticle contin- 
uum will be hard to detect experimentally. Averaging 
over the neutron spin polarizations, we obtain the dy- 
namic spin susceptibility x"(k, uj) as sum of 5 peaks with 
weights determined by various matrix element terms. 

In order to compare with the experiment of Ref. 
we add the contributions from horizontal and vertical 
charge modulations (Figs. Q and [21), and plot the result 
as function of the external momentum at fixed energy, 
furthermore we broaden the 5 peaks to account for the 
experimental resolution. 

In Figs. Hand [3 we show the response for bond- and 
site-centered stripe structures, where the parameters are 




55 meV 



36 meV 
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160 meV 



120 meV 




(7l,27t) 6 meV (271,71) 




(0,71) 



(71,0) 



FIG. 1: (color) Neutron scattering intensity, x"(k,a;), for a 
state with bond-centered charge order. Parameter values are 
c = 160 meV, Ai|(^^| = 100 meV, A2|<^:r| = 600 meV, A3 = 
A4 = 0. a = 550 meV is choosen is place the spin sector at its 
critical point. The panels show fixed energy cuts as function 
of momentum in the magnetic Brillouin zone. The 5 peaks 
have been replaced by Lorentzians with width F = 15meV , 
and the responses of horizontal and vertical stripes have been 
added. The figure can be directly compared to Fig. 2 of 
Ref.IH 



chosen to match the experimental result of Ref. 0. The 
coupling between charge and spin sector is captured by 
Ai which modulates the spin density, and in A2.3 which 
induce a spat ial variation in the exchange couplings. As 
in Refs. llalla . we observe a "dual" character of the lowest 
spin excitation branch: For small energies {e.g. 6 meV) 
the response consists of four peaks representing the four 
cones of spin-wave modes. With increasing energy the 
cones widen; however, the outer part becomes suppressed 
in intensity due to matrix element effects. Around 30 - 50 
meV the spectrum is dominated by the strong response 
near (7r,7r) ("resonance peak"), which arises from a sad- 
dle point of the mode dispersion. For higher energies, 
the modes gradually change their character towards a 
one-dimensional excitation spectrum, and the scattering 
intensity forms a diamond which moves outward with in- 



80 meV 200 meV 




(0,71) (7t,0) 



FIG. 2: (color) As Fig. but now for a state with site- 
centered charge order, and parameter values are Ai|<^a;| = 200 
meV, Aslt/i^l = 100 meV, A2 = A4 = 0, and s = 140 meV. 
Again the responses of horizontally and vertically ordered 
states have been added. 



creasing energy. Overall, there is reasonable agreement 
with the experimental data of Tranquada et al.^^ . A 
crucial point is the presence of two energy scales in the 
dispersion: a bandwidth of about 250 meV arising from 
the strong coupling along the stripes, and a saddle point 
at about 40 meV whose energy is dictated by the cou- 
pling across the stripes and the deviation of the order- 
ing wavevector from (7r,7r). Remarkably, the difference 
between the bond- and site-centered situation are mini- 
mal, i.e., the is little distinction on symmetry grounds. 
(Microscopically, however, strong spin fluctuations which 
drive the system away from the quasiclassical limit are 
favored in a bond-centered geometr jii^*iSi2S42i. ) 

Fig. 13 shows the result for bond-centered checkerboard 
order, with microscopic parameters similar to Fig. ^ 
Clearly, the result is completely different: The low-enery 
modes appear at four points rotated by 45 degrees com- 
pared to Figs. n El and the dispersion shows only a single 
energy scale. Thus, the "dual" character of the spectrum 
is absent here. 
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80 meV 200 meV 




(0,31) (3t,0) 



FIG. 3: (color) As Fig. ^ here for a state with two- 
dimensional plaquette charge order. Parameters are Ai|(/)a;| = 
100 meV, X2\(j}x\ = 600 meV, A3 = A4 = 0, and s = 1.3 eV. 
The spin order condenses at wavevectors (7r±7r/4, 7r±7r/4) in- 
stead of (7r±7r/4, tt), (tt, 7r±7r/4). Furthermore the dispersion 
features a single energy scale only, i.e., there is no magnetic 
response above the saddle point at (tt, tt). (A higher band of 
spin excitations start around 250 meV.) 



IV. CONCLUSIONS 

We have presented a general quantum lattice model 
which describes spin excitations in the presence of either 
static or fluctuating charge/bond order. 

In the case of anisotropic ('one-dimensional') charge 
order, the magnetic modes resemble semi-classical spin 
waves at low energies, but cross over into triplon exci- 
tations of a quasi-one-dimensional quantum paramagnet 
at higher energies. The crossover energy, associated with 
a saddle point in the mode dispersion, is only a fraction 
of the bandwidth for systems close to a magnetic quan- 
tum phase transition. Two ingredients are crucial for the 
"dual" character of the modes: (i) strong quantum fluc- 
tuations, and (ii) the presence of two energy scales in the 
dispersion. We found that the spectrum was relatively 
insensitive to the bond- or site-centered nature of the 
charge order. We also examined fully two-dimensional 
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'plaquette' ordered states, and found that they could not 
describe the experimental observations. 

An interesting open question is the precise link be- 
tween the neutron scattering observations and the mod- 
ulations observed in recent STM experiments. 
In principle, our approach can be adapted to any specific 
charge order observed in STM, and can then compute 
its spin excitation spectrum. It is already clear from our 
results that at least in La2-xBaxCu04, the charge or- 
der cannot be strictly a two-dimensional 'checkerboard' 
structure: anisotropic "quasi-onc-dimcnsional" correla- 
tions over some finite range appear to be required. 
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